Nucleation of superconductivity in mesoscopic star-shaped superconductors 
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We study the phase transition of a star-shaped superconductor, which covers smoothly the range 
from zero to two dimensions with respect to the superconducting coherence length £(T). Detailed 
measurements and numerical calculations show that the nucleation of superconductivity in this de- 
vice is very inhomogeneous, resulting in rich structure in the superconducting transition as a function 
of temperature and magnetic field. The superconducting order parameter is strongly enhanced and 
mostly robust in regions close to multiple boundaries. 
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PACS numbers: 73. 23.-b,74.25.Fy, 74.50. +r,74.80.Fp 



Many characteristics of a superconducting sample are 
strongly influenced by its size with respect to fun- 
damental superconducting length scales. Experimen- 
tal and theoretical investigations over the past few 
decades have clearly delineated the differences in be- 
havior between one-dimensional, two-dimensional and 
"bulk" three-dimensional samples in properties such as 
the critical field, critical temperature and critical cur- 
rent. Recently^, there has been renewed interest in the 
effect of sample dimension on the properties of supercon- 
ductors, fueled in part by their potential use in future 
nanometer scale devices. The complex geometry of such 
devices implies that one needs to consider the nucleation 
of superconductivity over a range of length scales in a 
single device. 

The nucleation of superconductivity in the presence 
of a magnetic field can be investigated by solving the 
Ginzburg-Landau(GL) equations for the superconduct- 
ing order parameter ip ancLthe vector potential A of the 
magnetic field B = V x AcT 
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Here a and (3 is the GL parameters, and n is the unit 
vector normal to the surface. 

Close to the critical temperature T c , |^| 2 is negligible, 
and one can analytically solve the linearized GL equa- 
tion ([!]) fojj some simple geometries. Saint- James and 
de Gennesa solved this linear equation for the case of 
a magnetic field parallel to the surface of a supercon- 
ductor. They found that superconductivity could exist 
at a field H c3 (the so-called surface critical field) larger 
than the upper critical field H C 2 in thepbulk supercon- 
ductor, H C 3 = 1.69.ff C 2. More recentlyETD, attention has 



focused on superconducting wedges subtending an angle 
7, with the magnetic field applied parallel to the wedge's 
edge. In the limit of small 7, it has been predictedEr 
that the surface critical field can be greatly enhanced, 
H C 3 = (1.73/7)if C 2, although this has never been demon- 
strated experimentally. This large enhancement of H c3 
can be ultimately traced to the increase in the "energy" 
(—a) in the GL equation ([!]) associated with the confine- 
ment of the order parameter. 

For more complicated geometries, and for regimes fur- 
ther away from the superconducting transition, it is nec- 
essary to solve the full GL equations numerically. This 
has been done for., many different .geometries, includ- 
ing wedges, ringstJll3, square loopsEaliJ and bridgesEj. 
Here we are interested in a four-pointed "star" geometry, 
shown schematically in Fig. 1. The remarkable feature 
of this geometry is that it encompasses a wide range of 
length scales. As we shall see, this results in appearance 
of regions of the sample with characteristically different 
behavior in the superconducting regime. 




FIG. 1. (a) Schematic representation of the method that 
was used to design samples, (b) Scanning electron micrograph 
of one of the samples. The inset shows magnified area of the 
apex. 

The experimental sample is prepared by removing por- 
tions of a circle of radius a from the corners of a square of 
side 2a, as shown in Fig. 1 (a). The experimental realiza- 
tion of this geometry is shown in Fig. 1(b). The sample 
is fabricated by conventional electron beam lithography 
with a 60 nm thick Al film. The apex-to-apex distance is 
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12 /im, the minimum dimension at each apex is less than 
100 nm. These dimensions are also used in all numerical 
calculations. In order to enable four terminal electrical 
measurements on the device, narrow electrical contacts 
of the same material are attached to opposite apices. 




FIG. 2. Evolution of \ip(x,y)\ 2 and magnetic field distri- 
bution in the star with temperature for an applied external 
field of H = 0.02# c (0) = 2 G, where H c (0) = 100 G is the 
bulk critical field at zero temperature, at different tempera- 
tures: a - T/T c = 0.99, b - 0.98, c - 0.97, d - 0.94. \ip{x,y)\ 2 
is normalized to its maximum value at these temperatures. 

To obtain the solutions of the GL equations (|l|), (||) 
with the boundary conditions (Bh, we use the finite- 



difference method applied earlier for the description of 
superconductivity in a mesoscopic square loopl3'Efl, us- 
ing the dimensions and parameters of the experimental 
sample of Fig. 1. One additional complication of this 
geometry for the finite-difference method is the infinite 
sharpness of the apices of the star, which makes the gen- 
eration of a mesh for the problem difficult. Fortunately, 
the vertices of the experimental sample have finite curva- 
ture, which results in a definite cut-off length of around 
100 nm for the tip of each apex, as shown in the small 
inset to Fig. 1. In our calculations of the order param- 
eter and magnetic field distributions, a square mesh has 
been used with the density of 800 nodes per side of the 
sample. At such a high density of nodes, the results of 
calculations are independent of the mesh. 

The solutions of the GL equations with appropriate 
boundary conditions permit us to obtain the spatial dis- 
tribution of the order parameter and the magnetic field as 
a function of temperature T, as shown in Fig. 2. As tem- 
perature is decreased slightly from T c (T/T c — 0.99), one 
observes nucleation of superconductivity as evidenced by 
a non-zero \if}(x, y)\ in two distinct regions of the sam- 
ple. The first region is a large area in the center, where 
one might intuitively expect superconductivity to nucle- 
ate. The second region is at the apices of the star, where 
superconductivity is enhanced by the close proximity of 
two boundaries, and resembles the situation that occurs 
in a superconducting "wedge." In between these two re- 
gions, |0| 2 is essentially 0, so that the superconducting 
phases of the sample are separated by regions of normal 
phase. The plot of the field distribution shows that the 
two superconducting regions are characterized by differ- 
ent behavior in a magnetic field: the central region shows 
the beginning of Meissner expulsion of the external field, 
while the magnetic field near the boundaries is very close 
to the external field value. As the temperature is low- 
ered further, the superconducting areas in the two regions 
grow as expected to cover almost the entire sample, and 
the normal regions between the apices and the center 
disappear. At the lowest temperature shown, the value 
of at the apices and in the central region achieves 
its maximum value corresponding to this temperature, 
\ip(x,y)\ 2 max - (1 - T/T c ) = 0.06. The magnetic field at 
the center is greatly reduced, showing the presence of a 
strong Meissner effect, but the magnetic field in a narrow 
region at the boundaries of the sample is still very close 
to the external field value. It should be noted that the 
value of \ip\ 2 at the apices is a maximum, even though 
the magnetic field is also a maximum. In fact, the value 
of | ■01 2 at the apices does not change appreciably over 
the temperature range shown in Fig. 2, demonstrating 
vividly the fact that superconductivity is the most ro- 
bust in regions close to multiple boundaries. 

The electrical contacts in the sample are placed at op- 
posite apices, so that the electrical current must traverse 
these apices, the central region, as well as any normal 
regions in between. Consequently, the inhomogeneous 
nature of superconductivity nucleation as the sample is 
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cooled through its transition imparts a distinctive shape 




FIG. 3. Experimentally measured resistive transition, 
taken at different magnetic fields. The inset shows the nu- 
merically calculated R(T). T c = 1.32 K, and R N = 12 CI 



to R{T), as shown in Fig. 3. The normal-to- 
superconducting transition starts with a decrease in resis- 
tance over a relatively broad range in temperature even 
at Hq = 0. In this temperature range, superconductiv- 
ity begins to nucleate at the apices and at the central 
part, as was shown in Fig. 2, but the regions near the 
apices, being narrower, contribute a large fraction to the 
resistance change. As the temperature is reduced fur- 
ther, there is a much sharper drop in the resistance when 
R is a little less than half the normal state resistance 
R 7y • Since the central region and the apices are already 
superconducting, this rapid change corresponds to the 
normal-to-superconducting transition of the "necks" be- 
tween the apices and the central region. Calculations of 
the superconducting transition confirm this picture. 

Numerically, we calculate the resistance of the sample 
using the order parameter distribution. As the order pa- 
rameter and magnetic field distributions, the resistance 
of the sample is calculated on a rectangular mesh. The re- 
sistance of a mesh cell in the sample is determined by the 
value of \ifj(x,y)\ 2 in the cell, being zero if \ifj(x,y)\ 2 > 0, 
and equal to the normal state resistance if \i/j(x, y)\ 2 — 0. 
One must also take into account the Josephson coupling 
between the superconducting regions, which reduces the 
resistance of the "necks". We suppose that a region of 
the normal metal of length < p ■ £ (where p ks 2 is an ad- 
justable parameter) between superconducting "islands" 
does not contribute to the resistance in the circuit. The 
inset of Fig. 3 shows result of this calculation for two 
different held values. The theoretical curves reproduce 
the main qualitative features of the experimental curves, 
including the "kinks" (denoted by the arrows) at lower 
resistance, and are in reasonably good quantitative agree- 
ment as well. The calculated position of the beginning 
of the resistive transition compares rather well with the 



experimental curves when the tunneling of the supercon- 
ducting electrons through the normal metal is taken into 
consideration. 

The behavior of the sample in a magnetic field (Fig. 
4) exhibits even more distinctly the difference between 
the nature of the superconductivity in the center and the 
apices. While \tp(x, y)\ 2 in the central region is rapidly 
attenuated with increasing magnetic field, its value in 
the apices does not show a large change. At large mag- 
netic fields, only the apices are superconducting. This 
would imply that the contribution to the resistance of 
the sample from the apices is not strongly affected by 
magnetic field. Evidence for this can be seen directly in 
the experimental transition curves. The resistance of the 
sample at the top of the transition, where the decrease in 
resistance is primarily due to the nucleation of supercon- 
ductivity in the apices, shows only a weak dependence on 
magnetic field. At lower values of resistance, where one 
has contributions from the central region as well as the 
"necks" , the field dependence is much stronger, reflecting 
the two-dimensional nature of the superconductivity in 
these regions. 

Even more striking is the magnetic phase diagram (T c 
vs. H) of the sample. Typically, this is obtained by 
varying the temperature to maintain the resistance of the 
sample at the midpoint of the transition (R/Rn = 0.5) 
while changing the magnetic field. Since our sample is 
inhomogeneous, and it is not clear how one would de- 
fine the unique T c , we have measured the phase diagram 
with the sample biased at various points R/Rn of the 
superconducting transition. Figure 5 shows the result 
of these measurements. Although the curves are differ- 
ent, each one shows two distinct regions, separated by a 
well-defined "kink" denoted by the arrow. At low fields, 
the curves show a quasi-linear dependence on Hq, which 
turns into a quadratic dependence at higher magnetic 
fields. It is well known that the critical temperature of 
a superconductor varies linearly with the magnetic field 
in two dimensions, and quadratically in one dimensionQ. 
The "kink" in each curve defines the crossover from two- 
dimensional behavior at low fields to one-dimensional be- 
havior at higher fields. This behavior is in agreement 
with the evolution of the superconducting regions of the 
sample shown in the numerical simulations of Fig. 3. 
At low fields, the central region and the "necks" of the 
sample are superconducting, giving a two-dimensional 
characteristic to the phase diagram. At higher fields, 
superconductivity in these regions is attenuated, leaving 
a finite order parameter only in the apices, which gives 
rise to a quadratic dependence on the magnetic field. It 
is interesting to note that the low-field linear behavior 
of all the curves extrapolates to a single temperature at 
Hq = 0, while the high-field quadratic behavior of all the 
curves also extrapolates to a single (but different) tem- 
perature. The different temperatures clearly point out 
the inhomogeneous nature of the superconducting tran- 
sition in this device. 
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FIG. 4. Evolution of \iji(x, y)\ 2 with magnetic field at T/T c = 0.8. \ip(x, y) 2 is normalized to its maximum value at this 
temperature. 
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FIG. 5. T — H phase diagrams measured at R/Rn = 0.58 
(a), 0.42 (b), 0.33 (c), 0.25 (d), 0.21^), 0.17 (f), 0.125 (g), 
0.08 (h). Inset shows the Little-Parkscj oscillations at smaller 
field for R/Rn = 0.04. 

In conclusion, we have analyzed the resistive transition 
in a star-shaped sample, which combines properties of 
two-dimensional and zcro-dimcnsional superconductors 
in a unique manner. Occurrence of qualitatively different 
regions in the resistive transition reflects a smooth change 
from 2D to 0D superconducting regime when increasing 
applied magnetic field or temperature. 
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